NASA TECHNICAL NOTE 


CO 

p««. 

NO 




LOAN COPY: RETURN a 
AFWL(DOUl) 
kirtland afb, n n m 


EVALUATION OF THE USE 
OF A SINGULARITY ELEMENT 
IN FINITE-ELEMENT ANALYSES 
OF CENTER-CRACKED PLATES 


by Alexander Mendelson, Bernard Gross, 
and John E. Srawley 

Lewis Research Center 
Cleveland, Ohio 44135 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, 



D. C. • MARCH 1972 

■ • ‘ 


TECH LIBRARY KAFB, NM 



tech library kafb, nm 

■■■linn 

0133407 


1. Report No. 2. Government Accession No. 

NASA TN D-6703 

3. Recipient's Catalog No. 

4. Title and Subtitle 

EVALUATION OF THE USE OF A SINGULARITY ELEMENT IN 
FINITE -ELEMENT ANALYSES OF CENTER-CRACKED PLATES 

5. Report Date 

March 1972 

6. Performing Organization Code 

7. Author(s) 

Alexander Mendelson, Bernard Gross, and John E. Srawley 

8. Performing Organization Report No. 

E-6680 

10. Work Unit No. 

134-03 

9. Performing Organization Name and Address 

Lewis Research Center 

National Aeronautics and Space Administration 
C leveland, Ohio 44135 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Note 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 

N5. Supplementary Notes 1 


16. Abstract 

Two different methods are applied to the analyses of finite width linear elastic plates with cen- 
tral cracks. Both methods give displacements as a primary part of the solution. One method 
makes use of Fourier transforms. The second method employs a coarse mesh of triangular 
second-order finite elements in conjunction with a single singularity, element subjected to ap- 
propriate additional constraints. The displacements obtained by these two methods are in 
very good agreement. These results suggest considerable potential for the use of a cracked 
element for related crack problems, particularly in connection with the extension to nonlinear 
material behavior. 


17. Key Words (Suggested by Authoi (s) ) 

Crack displacements 
Singularity elements 
Finite elements 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

23 

$3. 00 


18. Distribution Statement 

Unclassified - unlimited 


’ For sale by the National Technical Information Service, Springfield, Virginia 22151 



EVALUATION OF THE USE OF A SINGULARITY ELEMENT IN FINITE- 
ELEMENT ANALYSES OF CENTER-CRACKED PLATES 


by Alexander Mendelson, Bernard Gross, and John E. Srawley 

Lewis Research Center 

SUMMARY 

Two different methods are applied to the analyses of finite width linear elastic 
plates with central cracks. Both methods give displacements as a primary part of the 
solution. One method makes use of Fourier transforms. The second method employs 
a coarse mesh of triangular second-order finite elements in conjunction with a single 
singularity element subjected to appropriate additional constraints. The displacements 
obtained by these two methods are in very good agreement. These results suggest con- 
siderable potential for the use of a cracked element for related crack problems, par- 
ticularly in connection with the extension to nonlinear material behavior. 


INTRODUCTION 

The purpose of this study was to compare results of two methods of numerical anal- 
ysis that model the deformation of cracks in two-dimensional bodies of linear elastic 
materials. The study is preliminary to the further development of a hybrid finite- 
element method for application to cracked bodies of nonlinear materials. The hybrid 
method employs a special crack-tip element that can be designed to satisfy a variety of 
conditions. It was considered necessary first to evaluate the accuracy of the hybrid 
method for linear material behavior by comparison with an inherently more accurate 
method that makes use of Fourier transforms and satisfies compatibility and equilibrium 
throughout the body. The two methods were applied to centercracked plates of finite 
width for which linear elastic stress intensity coefficients have been obtained with high 
accuracy by several investigators (refs. 1 to 4). 

The Fourier transform method reduces the problem to the solution of a pair of 
Fredholm integral equations that can be solved by a simple iterative process which is 
rapidly convergent. The advantages of the method are that the crack surface displace- 


ments are a direct result of the solution of the integral equations, the stress (and strain) 
distributions are given by simple quadratures, and the required computational program 
is fairly simple. The ordinary finite -element method is inherently unsatisfactory for 
resolution of displacements very near to the crack tip because of its essentially discrete 
nature. The introduction of a finite element that itself contains a crack (ref. 5), called 
a singularity element (ref. 6), avoids this difficulty, but at the cost of loss of conditions 
of interelement continuity of displacement. In the present study a crack tip element of 
the inverse square root singularity type was incorporated into a coarse mesh of second- 
order elements to model the deformation of the crack and compute stress intensity coef- 
ficients. The results compare quite favorably with those of the Fourier transform 
method, and thus lend confidence in the further development of the hybrid method for ap- 
plication to nonlinear material behavior (which will require a different kind of crack tip 
element, however). 


SYMBOLS 


2a 

crack length 

B 

plate thickness 

2b 

plate width 

P>] 

elastic matrix for plane strain conditions 

FT 

Fourier transform method 

FE 

finite element method 

G 

shear modulus 

K I 

stress intensity factor, mode I 

M 

element stiffness matrix 

2L 

plate length 

N 

mesh density, number of finite elements per unit area 

{P 1 

nodal force column vector 

r, 0 

polar coordinates 

SyJ S x , S X y 

dimensionless stress 

[T] 

displacement matrix 

{U } 

displacement column vector 

u,v 

displacements in x and y directions 
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v o 

x,y 

[r] 

{e} 

e y ,e x’ ^xy 
X 

M 

v, I 


{crj 

v c 
{ } 

[ ] 


uniform end displacement of plate 

dimensionless cartesian coordinates with respect to b 

strain matrix 

strain column vector 

components of strain column vector 

plate length to plate width ratio 

Poissons ratio 

dummy integration variables 

stress column vector 

components of stress column vector 

column vector, m x 1 

matrix, m X m 


FOURIER TRANSFORM SOLUTION 

We consider the plane stress solution for an infinite strip loaded in tension. The 
strip is two units wide and contains a central crack of length 2a. Because of the sym- 
metry we need consider only one quadrant such as the upper right quadrant. The details 
of the solution are given in appendix A. Here, we summarize the equations and method 
used in the computations. If v(x, 0) is the displacement along the crack, then, as shown 
in appendix A, the following equations must be satisfied: 

/'a/b Z'a/b 

v(x, 0)=^ v(v,0)K di; + ^ [1 + dy (1) 

here, v is the displacement of the crack divided by the ratio of the load to the elastic 
modulas and 


where 



i^OO 

QU) = Jq Q(y)cos(|,y)dy 


( 2 ) 


( 3 ) 


3 


and 


pW - £ 


Q(|)K 2 (i/, $)d$ 


(4) 


The functions K, K^, Kg, and Ij are given by equations (19), (18), (25), and (21) in 
appendix A. 

The solution is obtained as follows. Assuming p( p) equal to zero, equation (1) is 
solved for v(x, 0). This is done by straightforward iteration since equation (1) is a 
simple Fredholm equation of the second kind, which converges very rapidly. Once 
v(x, 0) is determined the integrations indicated in equations (2) to (4) are carried out suc- 
cessively to determine p( p), which is then substituted into equation (1), and the process 
is repeated until there is no longer any significant change in either v(x, 0) or p(y). The 
whole process converges very rapidly, four or five iterations being sufficient for any 
engineering application. 

Once the crack surface displacement is determined from the solution of equations 
(1) to (4), the stress at any point can be computed by simple quadrature as follows: 


S x (x, y) = 1 - n 


/* a/b / 9i i \ r 


QU)K 3 d£ 


/■ aA ( ^ r 

S y (x, y) = -n I v{v, 0) flj + y — Jdp + J 


QUJKg d| 


V ( 5 ) 


/ a/b /•«> _ 

v(j/,0)I 2 di > + Q(?)K 4 d§ 


where the functions 1^, I 2 , K 2 , Kg, and K 4 are given in the appendix. We note that in 
performing the numerical integrations in equation (1), a difficulty arises from the fact 
that the kernel K(x, p) is singular at p = a/b. However, the principal value can be cal- 
culated as is shown in reference 7. 

Figure 1 shows results obtained by this method (originally published in ref. 4). On 
the scale of this plot the individual values of the dimensionless stress intensity coeffi- 
cient 2KjBb/Pa*/2 would be indistinguishable from the curve shown. Values taken 
from references 1 and 2 for comparison would also be indistinguishable from the curve. 
Further, results by this method are compared subsequently with results obtained by the 
finite -element method. 
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FINITE-ELEMENT ANALYSIS 


Consider the plane elasto- static solution for a plate of finite dimensions, subject to 
a uniform end displacement. The plate has a central crack of length 2a and width 2b 
(fig. 2). Because of symmetry, we need consider only one quadrant. A mesh is gener- 
ated of the type shown in figure 2(c). The stiffness matrices for the triangular and 
singular elements (fig. 3) are evaluated as shown in appendix B by the following equa- 
tions: 


[^Triangular = B // [r] T [D][r]dA (6) 

element 

Wsingular = B [ c ‘ 1 f \j[f [r^M^Ac' 1 ] < 7 > 

element 

The matrices [r], [D], [C], and [rj are defined by equations (28), (29), (36), and (37) 
in appendix B. 

The master stiffness matrix is then obtained by superimposing the various influence 
terms of the finite elements sharing common nodal points. The solution to the master 
stiffness matrix is then obtained by satisfying nodal equilibrium and boundary displace- 
ment conditions. 


NUMERICAL RESULTS 

The finite-element representation of the problem solved is shown in figures 2 and 3. 
The computational scheme is illustrated by the simple flow diagram as shown in fig- 
ure 4. A comparison of dimensionless displacements along the crack quadrant surface 
is given in figure 5. As shown in these figures, the finite -element displacements com- 
pare favorably with that computed with the more accurate Fourier transform method. 
Because the Fourier transform method satisfied the conditions of compatibility and 
equilibrium throughout the strip, one would expect that the results of this approach to 
be more accurate. It is noteworthy that the Fourier transform results correspond to 
elliptical crack profiles in all three cases. Table I lists the Y stress intensity coeffi- 
cients obtained herein and from references 3 and 5. Since displacements are the pri- 
mary solution, the Y values were based on the computed nodal displacements. As the 
ratio X decreased for a fixed number of elements (fig. 2), the comparison improved. 
This was primarily due to the fact that, as X decreases the finite elements per unit 
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area, N (or mesh density) increases. The Y stress intensity coefficients showed a 
much greater sensitivity to the X variation than the displacements. 

As shown in figure 3, the location of the crack tip was always taken at the midpoints, 
between nodes 3 and 5. Node 4 was always taken at the midpoint between the crack tip 
and node 3. A preliminary study revealed that the results obtained by locating the crack 
tip ±20 percent from the midpoint produced a small variation in the results. For larger 
changes in the crack tip location significant differences appeared. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 1, 1971, 

134-03. 
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APPENDIX A 


FOURIER TRANSFORM ANALYSIS 


We consider the problem of an infinite strip of width 2b and thickness B, contain- 
ing a central crack of length 2a. The strip is loaded in tension with a tensile stress 
cr^, normal to the plane of the crack. The stresses and geometric variables are non- 
dimensionalized in terms of a ^ and W, respectively. Because of symmetry, only 
one quadrant of the plate need be considered as shown in figure 6. The solution is ob- 
tained by the superposition of two solutions designated by A and B. 

Problem A consists of a plate with a central crack, a tensile load at infinity, a 
normal load p(x) along the crack surface, and a self-equilibrating side load Q(y) at 
x = 1 (fig. 6). Problem B consists of a plate without a crack and without a load at 
infinity, but with a side load -Q(y) at x = 1 (fig. 6). By superposing the solutions to 
both these problems with the proper choice of Q(y) one obtains the desired solution to 
the problem of figure 6. 

The solution to problem A can be obtained as follows. The equilibrium and com- 
patibility equations for the plane problem of elasticity are 


3S 9S 

_JS + _22 = o 

3x 9y 


!!z + !^ = 0 

3y 9x 




( 8 ) 





Applying a Fourier sine transform to the first of these equations and Fourier cosine 
transforms to the last two, solving the resulting ordinary differential equations and then 
taking the inverse transforms results, as shown in detail in reference 7, in 
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S v = 1 + 2 ^ A m (l + ni7ry)e" m7ry cos mTrx 




m=l 


S = 2 ^ A (1 - m77y)e~ m7ry cos m7rx 

x m=l m 




OO 

s = 27ry mA e _m7ry sin m 7 rx 

xy m=l m 


J 


( 9 ) 


That equations (9) indeed satisfy equations (8) can be verified by direct substitution. 

The displacements can now be calculated from the usual stress-strain relations and 
strain displacement results resulting for the case of plane stress in 


v = 



v(i/,y)d*/ - - 

7T 



A 

— — [2 + (1 + /i)mjry]e _m7ry cos mTrx 
m 


( 10 ) 


Along the plane of the crack (y = 0) 


and 



(ID 


Since equations (11) represent Fourier series expansions, it follows that 


and also 


A 


m 



SyU, 0)cos mrr^ d| 


m 



v(i>, 0)cos mirv dy 


( 12 ) 


(13) 
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The upper limit of integration in equation (13) follows from the fact that 


v( y , 0) = 0, 


- < x < 1 


Assuming a normal stress distribution p(x) along the crack surface, equation (12) can 
be written 

/ a/b /*1 

p(£)cos mffS d£ + J o (|,0)cos mnE, d£ (14 

v *Vb y 

Now, substituting (14) into the second of equations (11) and equation (13) into the first of 
equations (11) results in 


/ a/b i /*a/b , 

\(v,0)d v + - / p(i/)ln| 

77 •'O 


(cos nv - cos 


nx)^] d i 


- f S (£, 0)ln [4(cos - cos 77x)^]d(; 
77 •'a/b y 


where 


S (|,0) = 1+1 / v(y,0) L- .. ? . 9 s . JT^ c0 - s dy 


(cos 77 v - COS 77^)' 


and where the following identities have been used 


00 

I: 


cos mnv cos m77X = - - In [4(cos nu - cos ttx)^! 

4 


£ m cos m77z, cos - * — cos EE £2g H 

m= l 2 ( cos nv _ cos 7^)2 


The first identity in equations (17) can be found for example in reference 8 (p. 358). 

The second can be obtained by differentiating the first one twice. It is obvious of course 


9 



that the second sum is divergent. Whether the formal sum given has meaning can only 
be ascertained from the final answer, which in this case can be shown by direct substitu- 
tion to satisfy the equations of equilibrium and compatibility. 

Finally substituting equation (16) into (15) gives 



2 2 


+ 



2 tt sin 


7r(x + y) 


2 2 


(19) 


For a given pressure distribution p(j/) along the crack, equation (1) represents a Fred- 
holm equation of the second kind. It can be solved by a straightforward iteration which 
converges very rapidly. 

Once the crack opening is determined by the solution of equation (1), the stress dis- 
tribution can be obtained by direct numerical integration as follows: Substituting equa- 
tion (13) into equation (9), reversing the order of integration and summation, and for- 
mally summing the resulting infinite series give 
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( 20 ) 


/' a/i) / 3IA 

S x (x, y) = -7T J V{v, 0) ^ + y — ijdy 


S xy (x, y) = Try 


/ a/b 


31, 

v(y, 0) — - d v 

dx 


where 


I = l j cos + x)cosh 7ry - 1 cos tt(^> + x)cosh 7ry - 1 

1 A I O < 


2 2 
[cos tt(i^ + x) - cosh 7iy] [cos n(v - x) - cosh ;ry] 


( 21 ) 


At the side boundary x = 1 we now have a self -equilibrating normal load 


Q(y) = s x ( i, y) 

which can be computed from the second of equations (20). We have thus obtained a so- 
lution to problem A. 

The solution to problem B is given in reference 9 (p. 412). For an infinite strip 
with a symmetric edge load Q(y), we take the infinite Fourier transform 


/*oO 

Q U) = Jq Q(y)cos(£, y)dy 


( 22 ) 


Then the stresses are given by 


_ 4 / (sinh | + £ cosh Qcosh gx - gx sinh g sinh gx CQg ^ y ^ 

JT J Q 2g + sinh 2g 


"N 


g _ 4 / (sinh g - $ cosh g)cosh gx + gx sinh g sinh gx CQg ^ y 

*«/) 


r (23) 


2g + sinh 2g 


xy 


s _ = i / q( ?) sinh 1 cosh gx - j cosh j sinh gx gin |y d| 

'Jo ’ 


2g + sinh 2g 




Now at y = 0 this solution gives a stress in the y direction given by the second of 
equations (23) 


P(y) = S (x, 0 ) = f 

y o 


g)dg 



where 


4 (sinh g - g cosh g)cosh ^y_+_^y_ sinh g sinh 
ti 2g + sinh 2g 


(24) 


(25) 


By combining the solutions of problems A and B we can obtain a solution in which both 
the side loads Q(y) and the normal loads on the crack surface cancel. This can be done 
by the simple iterative procedure described in the body of the report. 
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APPENDIX B 


FINITE ELEMENT ANALYSIS 


In the following analysis it is assumed that the body forces are negligible and that 
the model material is homogeneous and isotropic, has a constant thickness, and is sub- 
jected to plane strain elastostatic conditions. Define 


{ U} 





U 1 

u 2 





Thus, we can also write 


{ U} = [T] { 5 } 

{e} = [r]{ 6} 

{a} = [D] { e } = [D][r]{ 6} 


(27) 

(28) 
(29) 


Imposing virtual displacements at the element nodal points and equating the work done by 
by nodal force { P} to the internal work done by the stress { a} (ref. 10), we obtain 

|6 v J T {P} =B j //'|e v } T {a}dA 
area 


where superscript v refers to virtual and T to transpose. Substituting equation (29) 
and (28) into the last equation and comparing coefficients yields 


13 



{P} = 



[r] [D][r]dAl {6} 


Hence, 


{P}=[k]{6} (30) 

where 

[k] = B If [r] T [D][r]dA (31) 

area 

The typical triangular element as shown in figure 3 is assumed to be a second-order 
element. The displacement functions are as follows: 

6 

v(x, y) = V W m (x, y)v m 
m=l 

6 

u(x, y) = V W m (x, y)u m 
m=l 


where 


W m^ x> y) = a l, m + a 2, m x + a 3, m y + a 4, m x + a 5, m 3 ^ + a 6, m y 2 


The constants a- i = 1, 2, . . . , 6 are evaluated such that, at location (x ,y ) of 
1 3 m 7 \ J txi 

node m, W m (x m ,y m ) = 1 and W m (x^,y^) =0 at the remaining l nodal points. Thus 
[T(x, y)] is obtained. The element stiffness matrix is evaluated in a manner similar to 
that in reference 11. 

For the singularity element as shown in figure 3 the displacement functions, taken 
from references 12 and 13, are assumed to have the general form 


v(r, B) 


u(r, 9 ) 


2G 

_ 1 _ 

2G 


’ 12 


E v n/ VM> 


n=l 

> 

*12 


E V n/2g n {n ’ 9) 



J 


14 



For odd values of n 


where 


and 


f n (n, 0) = Fj(n, 0)sin 9 + F 2 (n, 0)cos 9 
g n (n, 9) = F^n, 0)cos 9 - Fg(n, 0)sin 9 


Fj(n, 9) = Q - n - 4/^ cos^n - ^0 + |n - ~^ cos ^ + 


F 2 (n, 9) - ^ + n - 4/jj sin^n - 1^0 - - -^sin^n + -^0 


for even values of n 


f n (n, 9) = Fg(n, 0)sin 0 + F^(n, 0)cos 0 
g n (n, 0) = F 3 (n, 0)cos 0 - F 4 (n, 0)sin 0 


(32) 


(33) 


(34) 


where 


and 


Fg(n, 0) = (3 - n - 4/i)cos(n - 1)0 + (n + l)cos(n + 1)0 


> 

F 4 (n, 0) = -(3 + n - 4^t)sin(n - 1)0 + (n + l)sin(n + 1)0 


(35) 


To obtain the singularity element stiffness matrix, the element is divided into many (p) 
segments. A typical segment is shown in figure 3. From equation (27) in matrix formu- 
lation we have 


{ U} = [H]{ d} 

At the six nodal points we evaluate 

{6} = [C]{ d} 


(36) 
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hence 


{d} = [C -1 ]{ 6} 


thus 

{U} = [H][C” 1 ]{ 6 } 

And from equation (28) 

{€} = [r 1 ]{ d} = [r 1 ][C" 1 ]{ 6} (37) 

and it follows that 

[^singularity = [ff Ir/pIr^CT 1 ] (38) 

element 


Let 


[P q t ] = r T ][D][r]dA 

4 ’ 1 nrpn 


After some algebraic manipulation and integration the matrix elements become 

j=l s=l z=l 4 


(39) 


where q, t = 1, 2, . . . , 12, p is the number of segments of the cracked element, and 
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G (®)l,2n-1 = F l( n > ®H x i sin2fl + cos 2 e) + F 2 (n, 0)X 2 sin 9 cos 9 + F'j(n,9)sin 9 cos 9 + F 2 (n, 9)cos 2 9 
g (^ 2, 2n-l = F l( n ’ 0 K x i c °s 2 0 + sin 2 9) - F 2 (n, 9)Xg sin 9 cos 9 - Fj(n, 9)sin 9 cos 9 + F 2 (n, 9)sin 2 9 
°( e h, 2n-l = F l( n ’ °'> X 2 sin 29 + F 2^ n > 9 > x 2 cos 29 + F i< n > e ) cos 29 ~ F 2 ,( n > 9 ) Bin 28 

0 o 

2n = F 3< n > 9 K n sin 9 + cos^9) + F 4 (n, 9)Xg sin 9 cos 9 + F^(n, 9) sin 9 cos 9 + F^n, 9)cos 2 9 
G ^2, 2n = F 3^ n ’ 0 ^ n cos2 ® + sin 2 9) - F 4 (n, 9)Xg sin 9 cos 9 - Fg(n, 9) sin 9 cos 9 + F 4 (n, 9)sin 2 9 
G ^ e ^3, 2n = F 3^ n > e ) X 3 sin 29 + F 4^ n > ®) X 3 cos 29 + Fg(n, 9)cos 29 - F 4 (n, 9)sin 29 

where 

X t = n - - 

1 2 

v _ _ 3 

1 2 

X 3 = (n - 1) 

— [F(n, 0)] = F’(n, 0) 
d0 


and 
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TABLE I. - DIMENSIONLESS STRESS INTENSITY 


COEFFICIENTS FOR VARIOUS NOTCH 
DEPTH TO PLATE WIDTH RATIOS 


AND MESH DENSITY 


Method 

Source 

Mesh 

density, 

N 

Notch depth to plate width 
ratio, a/b 

0. 1 

0.3 

0.333 

0.305 

Dimensionless stress in- 
tensity coefficient, 

Y = (2bBK I )/(Pyi) 

Fourier 

This report 

-- 

1. 78' 

1.87 

1.90 

2. 11 

transform 

Ref. 3 

-- 

1.79 

1. 88 

— 

2. 12 

Finite ele- 

Ref. 5 

(a) 

— 

— 

1.83 

— 

ment 

This report 

16 

— 

— 

1.74 

— 




18 

— 

— 

— 

2.08 




21 

— 

1.79 

1. 81 

— 




32 

1. 52 

1.82 

— 

— 




64 

1.66 

— 

— 

— 




80 

1.65 

— 

— 

— 


a Not available. 



Figure 1. - Stress calibrations for center-cracked specimen 
(refs. 1, 2, and 7). 


CD 



v 0 





(a/b = 0.100, or 0.300, 0.333) 



(a) Center cracked plate subjected (b) Dimensionless (c) Typical finite-element mesh for quarter section, 

to uniform end displacement, Vq. quarter section. 

Figure 2. - Plate configuration and typical finite-element mesh. 
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Crack element with node locations at Typical 2nd order triangular element with node 

verticies and midpoints and typical locations at verticies and midpoints 

segment used for numerical inte- 
gration 


Figure 3. - Typical finite-elements for center cracked plate. 



Figure 4. - Simplified flow diagram of program operations. 
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(a) Notch depth to plate width ratio, 0. 1. 
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Dimensionless location along crack, -x/(a/b) 


(b) Notch depth to plate width ratio, 0. 3. 

Figure 5. - Crack displacement profile for various length to width ratios. 



Dimensionless location along crack, -x/(a/b) 

(c) Notch depth to plate width ratio, 0. 5. Length to width ratio, 2. 
Figure 5. - Concluded. 



Figure 6. - Fourier transform solution by superimposing problem A and problem B. 
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